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ABSTRACT 

If X-ray flashes are due to the forward jet emissions from gamma ray bursts 
(GRBs) observed with large viewing angles, we show that a prompt emission from 
a counter jet should be observed as a delayed flash in the UV or optical band 
several hours to a day after the X-ray flash. Ultraviolet and Optical Telescope 
on Swift can observe the delayed flashes within ~ 13 Mpc, so that (double-sided) 
jets of GRBs may be directly confirmed. Since the event rate of delayed flashes 
detected by Swift may be as small as ~ 6 x 10~ 5 events yr _1 , we require more 
sensitive detectors in future experiments. 



Subject headings: gamma rays: bursts — gamma rays: theory 
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1. INTRODUCTION 

Several observations suggest that gamma-ray bursts (GRBs) are caused by relativistic 
jets (e.g., Frail et al. 2001). However, in order to establish that GRBs are collimated, 
other observations are indispensable, such as polarization observations (Ghisellini & Lazzati 
1999; Sari 1999) and microlensing observations (Ioka & Nakamura 2001b). Some theoretical 
models of jet emissions have been discussed (Totani & Panaitescu 2002; Huang, Dai & Lu 
2002; Dado, Dar & De Rujula 2001). If GRBs are due to forward jet emissions, there should 
most likely be counter jet emissions, as in the AGN (Begelman, Blandford & Rees 1984) 
and the microquasar (Mirabel & Rodriguez 1999). Therefore the detection of counter jet 
emissions will give us direct evidence for the jet model of GRBs. 

The confirmation of a counter jet has been by far the most important factor in the 
jet model of astrophysical objects. A mysterious spot was found in SN1987A using the 
speckle technique (Meikle et al. 1987; Nisenson et al. 1987). Many models including the jet 
model were proposed (Rees 1987; Piran & Nakamura 1987). At that time, it was difficult to 
distinguish each model from observations since only one spot was found. In the jet model, 
the counter jet should be observed although it is dim due to redshifting (Piran & Nakamura 
1987). However, later in 1999, two spots were confirmed using new software to analyze the 
speckle data (Nisenson & Papaliolios 1999). Very recently the jet feature of the ejecta of 
SN1987A whose position angle is the same as the mysterious spot was confirmed by the 
HST (Wang et al. 2002). As a result the jet model by Piran & Nakamura (1987) took 
the advantage. Furthermore the observation of a counter jet may enable us to estimate the 
Lorentz factor of the jet, as for the AGN and microquasar. Therefore it is important to 
argue the observational properties of the emission from the counter jet of a GRB. 

Let us consider the emission from a counter jet with a Lorentz factor 7. The observed 
typical frequency of the counter jet emission is about 7 2 times smaller than that of the 
forward jet (i.e., the GRB). Since the typical frequency of the GRB is ~ 100 keV, the typical 
observed frequency of the counter jet emission is ~ 10(7/100)~ 2 eV, which is in the UV or 
optical band. This transient phenomenon should be observed about several tens of hours 
after the forward jet emission, since it is at a radius of order 10 14 -10 15 cm that photons are 
emitted from each jet leaving, almost simultaneously, the central engine. We call this event 
the delayed flash (DF). 

Any attempt to detect the DF might be difficult since the afterglow of the forward jet 
might be brighter than the DF. The afterglow of the GRB, i.e., the afterglow of the on-axis 
forward jet, is much brighter than the DF. However if the forward jet is observed with a 
large viewing angle, there is a chance to observe the DF since the forward jet emission is 
also dim at the time of the DF. 
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Recently, we studied the emission from the off-axis jet (Yamazaki, Ioka, & Nakamura 
2002a, b; see also Nakamura 2000; Ioka & Nakamura 2001a). We proposed that if we observe 
a GRB with a large viewing angle, it looks like an X-ray flash (XRF), a new class of X-ray 
transients which has been recently recognized as a phenomenon related to the GRB (e.g., 
Heise et al. 2001; Kippen et al. 2002; Barraud et al. 2002). The off-axis jet model can 
explain the typical observed frequency and other observational characteristics of the XRF, 
such as the peak flux ratio and the ffuence ratio between the 7-ray and the X-ray band, the 
X-ray photon index, the typical duration, and the event rate (Yamazaki, Ioka, & Nakamura 
2002a, b). Although the origin of XRFs is uncertain, we assume that XRFs arise from 
prompt off-axis jet emissions hereafter. 

In this paper, we will show that the DF can be observed after the XRF in principle. We 
will calculate the light curves of the XRF, DF, and the afterglow of the XRF, and discuss 
whether the DF can be detected by the Swift satellite. We will find that we need more 
sensitive detectors to detect the DF. In § 2, we describe a simple forward-/counter-jet model 
for the XRF and DF. In § 3 and § 4, we show the light curves of the XRF, the DF, and the 
XRF afterglow. § 5 is devoted to a discussion. 

2. INSTANTANEOUS EMISSION FROM AN EXPANDING JET 

We extend the simple jet model by Ioka & Nakamura (2001) and Yamazaki, Ioka, & 
Nakamura (2002a). In these works we assume that the shell width / is much smaller than the 
separation of the shells L in the internal shock, since the separation L mainly determines the 
emission timescale of the forward jet. However in this paper we consider a finite shell width, 
since the shell width I determines the emission timescale of the counter jet. The cooling 
timescale is much shorter than other timescales (Sari, Narayan & Piran 1996), so we assume 
an instantaneous emission at the shock front as before. We use a spherical coordinate system 
(t, r, 9, 0) in the Lab frame, where the 9 = axis points toward the detector, and the central 
engine is located at r = 0. The forward jet has a viewing angle, 9 V > 0, which the axis of the 
emission cone makes with the 9 = axis, while the counter jet has a viewing angle 9 V + n. 
When the emitting shock front moves radially from t — to and r = tq with the Lorentz factor 
7 = 1/ yjl — (3 2 , the emissivity for the XRF has a functional form of 



where f(v') represents the spectral shape. The Heaviside step function H(x) describes that 
the emission is inside a cone of an opening half- angle A9. Then the observed flux per unit 
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frequency of a single pulse at the observed time T is given by 

F(T) 2r 2 7 2 f h 2 (l-(3 cos 0(T))] A0(f )/[z/ 7 (l - g cos 0(f))] 

^ J (3D*(r /c(3) J W [ 7 2 (l-/3 cos 0(f))] [ 7 2 (1 - /3cos0(f))] 2 ' lJ 

where 1-/5 cos 0(T) = (cP/r )(T-T ), l-/3cos0(f) = [1- (3 cos 9 (T)]/[(c(3/r )(t-T )\ and 
T = f -r /c/3 (see Ioka & Nakamura 2001a). For the XRF, 0(f) varies from max{0, 0„- A0} 
to 0„ + A0. The function A0(f) is given as 

/ 7T (0„ < A0 and < 0(f) < A6 - 9 V ) 

A ^) = \ rcosA^cosMcos^] (otherwige) • (3) 



COS 



sin 9 V sinS(t) 



Also for the DF we can use equation (2). For the DF, 0(f) varies from ix + max{0, 9 V — A0} 
to 7r + 9 V + A0, and the function A0(f) is given as 

J 7T (0„ < A0 and 7T < 0(f) < 7T + A0 - 0„) 

A ^) = i r EgiAg±£2^)^ l (otherwise) ' (4) 



cos 



— sin 9 V sin 9(t) 



The normalization of emissivity A(t) is determined by the hydrodynamics. Here for 
simplicity we adopt the following functional form, 

A(t)=A ( t -^%) 2 H(t-t Q )H(t e -t), (5) 
V r o/ c P/ 

where the emission ends at f = f e and the released energy at each distance r is constant. 
Our conclusion does not depend on f e or the functional form so much. A pulse-starting time 
and ending time are given as 

T^rP = To + (r o /c/3)(l-/3cos(max[O,0,-A0])), (6) 

T^ d RF) = r o + [(r o /c/3)+f e -fo](l-/3cos(0, + A0)), (7) 

for the XRF, and 

TL D a rt = T o + (r o /c/3)(l + /3cos(0, + A0)), (8) 

TenP = T o + [(ro/c/3)+f e -fo](l + /3cos(max[O,0,-A0])), (9) 

for the DF. 

The spectrum of the GRB is well approximated by the Band spectrum (Band et al. 
1993). In order to have a spectral shape similar to the Band spectrum, we adopt the 
following form of the spectrum in the comoving frame, 

f(u') = (u'/u' ) 1+ ^(l + u'/u' o r—, (10) 
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where cub (Pb) is the low (high) energy power law index. In the GRB, as ~ — 1 and (5b ~ —3 
are typical values (Preece et al. 2000). Equations (2), (5) and (10) are the basic equations to 
calculate the flux of a single pulse, which depends on following parameters: 7 > 1, 9„ < 1, 
A9 < 1, ji/ Q , r /c(3-f 2 , a B , Pb, D, A , t and t e . 

Hereafter we consider mainly the following canonical set of parameters; 7 = 100, r yA9 = 
5, r //9c7 2 = Is, as = —1, Pb = _ 3, /t/z/q = 200 keV and t = r /cp. We adopt t e = 
1.3r /c/3 since most pulses rise more quickly than they decay (Norris et al. 1996). The value 
/ yA9 = 5 has been obtained from the fitting of the afterglow light curve (Frail et al. 2001; 
see also Panaitescu & Kumar 2002). We fix the amplitude Aq so that the isotropic 7-ray 
energy E iso = AttD 2 S(20 - 2000 keV) satisfies (A9) 2 E iso = 1 x 10 51 ergs when j9 v = (Frail 
et al. 2001). Here S(vi - v 2 ) = J T e " d t F(T; v x - u 2 )dT is the observed fluence in the energy 
range v x — v 2 and F(T; v x — v 2 ) = f" 2 F u (T)du is the observed flux in the same energy range. 
Then, we obtain A = 0.24 erg s~ 1 cm -2 Hz -1 for the fiducial parameters. Note that the 
observed flux is proportional to D~ 2 } 



3. LIGHT CURVES OF X-RAY FLASH AND DELAYED FLASH 

In this section, we plot the light curves of the XRF and the DF using Eq.(2), and discuss 
whether these events can be observed by Swift satellite. 

First, we show the light curves of the XRF F(T; 15 — 150 keV) in Figure 1 with varying 
r )9 v . The observation band corresponds to that of the Burst Alert Telescope (BAT) on Swift. 
As 70„ increases, the peak flux of the XRF F^ a f F ^ decreases due to the relativistic beaming 
effect. 

The light curves of the DF F(T; 1.9 — 7.3 eV) are shown in Figure 2 with varying r )Q v . 
The observation band corresponds that of the Ultraviolet and Optical Telescope (UVOT) 
on Swift. We find that the flux remains almost constant in each pulse, and that the peak 
flux F^ff) does not depend on the viewing angle j9 v so much. This is because the value 
of 9(t) ranges between ir + max{0, 9 V — A#} and n + 9 V + A9 and (c/3/r )(t — T ) ~ 1 so 
that 9(t) ~ 9{T) ~ 7r in equation (2). From equation (2) the peak flux of the DF can be 
estimated as f££ ~ AvF v ~ Av(2r 2 -f 2 / (3D 2 )A A(j)f /(2-f 2 ) 2 ~ 10" 19 ergs s" 1 cm" 2 , where 
Au ~ 10 14 Hz, A0 ~ A9/9 V ~ 0.1 and / ~ 0.2. 



1 When we consider the effect of cosmology (Qm = 0.3, Cl\ = 0.7, and h = 0.7), D ~ 1 Gpc corresponds 
toz~ 0.2. Since we consider the case D < 1 Gpc in the following sections, this does not affect our argument 
qualitatively but alters the quantitative results up to a factor of 2. 
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The limiting sensitivity of the UVOT (BAT) can be estimated as 1 x 10~ 15 ergs s _1 cm~ 2 
(5 x 1CT 10 ergs s -1 cm -2 ) for a duration of ~ 5 x 10 3 (~ 10 2 ) seconds. The BAT localizes the 
XRF and the following observation by the UVOT may identify the associated DF. One can 
see that the DF with D < 13 Mpc is observable. Then the BAT can detect the preceding 



The start and end time of the DF is about Tg tart ~ 2t ~ 2 x 10 4 s and T end ~ 2t e ~ 
2.6 x 10 4 s for 7 = 10 2 , to — r o/ c /3 = 10 4 s and t e = IMq. Therefore, one should compare 
the flux of the DF with that of the XRF afterglow. In this section, we plot the light curves 
of the XRF afterglow and see whether or not the DF can be detected. We use model 1 of 
Granot et al. (2002) as a simple model of the off-axis afterglow emission from the collimated 
jet (see also Dalai, Griest, & Pruet 2002). 

For V = 0, the standard afterglow model, i.e., the synchrotron-shock model, can explain 
observational properties of the GRB afterglow very well (Piran, 1999), and gives the observed 
flux per unit frequency as F U (T; V = 0) = i^ R ~ SPH ) = G(u,T), where i^ R_SPH ) [ s the 
observed flux given by Rhoads (1999) and Sari, Piran & Halpern (1999). For V > A0, the 
emission is assumed to be from a point source moving along the jet axis. Then the flux is 
given by F V (T\ V ) = cPGfa-^aT), where a = (1 - (3)/(l - f3 cos0) ~ (1 + (70) 2 ) -1 . We 
choose = max(0, V — A0) to make this simple model more realistic (Granot et al. 2002). 
The Lorentz factor of the shell 7 can be determined by 



where tj et = 1.9 x 10 4 secn~ 1 / 3 (A^/0.05) 8 / 3 (£'/2 x lO^ergs) 1 / 3 is the jet-break time observed 
from an on- axis observer 2 , where E is the isotropic equivalent value of the total energy in 
the shock. We assume E = ri~ l E iso with a constant factor 7] 7 = 0.2, which is adopted in 
Frail et al. (2001). Then, we obtain E = 2 x 10 54 ergs and the geometry corrected total 
energy in the shock (A0) 2 E/2 = 2.5 x 10 51 ergs. 

Using above equations, we calculate the light curves of the afterglow. In order to study 
the dependence on the viewing angle V , we fix the rest of the parameters: the power- 
law index of accelerated electrons p = 2.25, the number density of the ambient matter 



XRF if 7^ < 30. 



4. 



AFTERGLOW OF X-RAY FLASH 




(11) 



2 For simplicity, we assume the relation t' = i?/(47 2 c), where t' and R are the time measured by an on-axis 
observer and the radius of the shock. 
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n = 1 cm -3 , e e = 0.1 and eb = 0.01, and the distance D = 1 Gpc. Figure 3 (Figure 4) shows 
the result in the case of A9 = 0.05 (A9 = 0.1). The observation band is 1.9-7.3 eV, which 
corresponds to that of UVOT. 

We also plot the UV flux of the DF in the same figures. We can see that for the 
canonical set of parameters (A9 = 0.05, 7 = 100, and r /c/37 2 = 1 sec), the UV flux of the 
DF dominates the afterglow when 9 V > 0.21. For comparison, we show the light curves of 
the DF with one of parameters changed from the fiducial value. For large 7, it is difficult to 
detect the DF since the starting and ending time of the DF is late and the flux of the DF 
is low due to the strong beaming effect. When we alter ro/c/3'-f 2 , the starting (and ending) 
time and the flux of the DF have a dependence oc (r Q /c/3'y 2 ) and oc (r /c/57 2 ) -1 , respectively. 
So the large r /c/37 2 case has qualitatively the same behavior as the large 7 case. One can 
easily find that when A6 becomes large, the flux of the afterglow of the XRF becomes large 
while the light curves of the DF remains almost unchanged. Therefore, we can conclude 
that the DF from a jet with smaller A6, 7 and ro/c/97 2 have larger chance to be seen. In 
consequence, according to the off-axis jet model, it is preferable for the detection of the DF 
that the preceding XRF has a low peak energy of a few keV, a small variability owing to 
large 9 V (Yamazaki, Ioka, & Nakamura 2002a,b) and a short duration due to small r /c/57 2 . 

We have used in this section model 1 of Granot et al. (2002). A more realistic model 
for the off-axis emission from the forward jet, model 3 of Granot et al. (2002), may have 
a more moderate rise before the peak of the observed light curve than the model we have 
adopted. However, we consider the case in which the viewing angle is as large as 9 V > 5A0, 
so that the differences between these models may be small. 

5. DISCUSSION 

We have calculated the light curves of the DF, XRF, and the afterglow of the XRF. 
We have shown that in principle, the DF emission can be seen in the UV band about 10 4 - 
10 5 seconds after the XRF if the viewing angle is large enough (about 0.2-0.3 rad) for the 
afterglow of the XRF to be dimmer than the DF. Since the UV flux of the GRB afterglow 
is much larger than that of the DF, only the DF associated with an off-axis jet, i.e., an 
XRF has any chance to be observed. The preceding XRF should have a low peak energy 
of a few keV, a small variability, and a short duration for the DF to be detected. Due 
to the relativistic beaming effect, the flux of the DF is so small that only nearby events 
(< 13 Mpc for the canonical parameters) can be observed by UVOT on Swift. Following 
Yamazaki, Ioka, & Nakamura (2002a), we can roughly estimate the event rate of the DF 
for the instruments on Swift as Rdf ~ 6 x 10~ 5 events yr -1 , where we adopt the event 
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rate of the GRBs r grb = 5 x events yr -1 galaxy -1 and the number density of galaxies 
n g = 10~ 2 galaxies Mpc~ 3 . Therefore, we need next-generation detectors, which are more 
sensitive than the instruments on Swift, to detect the DF associated with very dim XRF 
more frequently. 

The DF may be obscured by dust extinction. In fact, about half of accurately localized 
GRBs do not produce a detectable optical afterglow (Fynbo et al. 2001; Lazzati, Covino, 
& Ghisellini 2002). One explanation for these "dark GRBs" is that most GRBs occur in 
giant molecular clouds (e.g., Reichart & Price 2002). In this picture, a GRB has a detectable 
optical afterglow only if the burst and the afterglow destroy the dust along the line of sight 
to the observer (Waxman & Draine 2000; Fruchter, Krolik, & Rhoads 2001), as suggested 
by the comparison between X-ray and optical extinction (Galama & Wijers 2001). In this 
case the DF is obscured since the flux of the XRF is too dim to carve out a path for the DF. 
However this picture may have some problems, such as no evidence of an ionized absorber 
(Piro et al. 2002) and variable column density in the X-ray afterglow (Djorgovski et al. 
2001c). There are other explanations for dark GRBs, such as high redshift effects, dust 
extinction in the interstellar medium of the host galaxy (Ramirez- Ruiz, Trentham, & Blain 
2002; Piro et al. 2002) and so on. Therefore at present we cannot conclude that the DF is 
obscured. 

If we assume that the absolute magnitude of the host galaxy is about ~ —20 mag 
(Djorgovski et al. 2001a, b), the apparent magnitude is about ~ 20 + 51ogDGpc- Since a 
host galaxy with a size ~ 10 kpc has an angular size of ~ 10-D,^, arcsec, we can observe a 
point source which is dimmer than the host galaxy by ~ 10~ 4 -Dq pc if the angular resolution 
is ~ 0.1 arcsec. Therefore the DF has to be brighter than ~ 30 mag, and we can observe 
the DF if D < 13 Mpc. 

If the GRB is associated with a supernova (SN), the emission from the SN may hide the 
DF. The UV flux of SN1998bw was about ~ 17 mag at the distance D ~ 40 Mpc (Galama 
et al. 1998), i.e., ~6x 10~ 15 D^c er S s_1 cm ~ 2 > so that a SN like SN1998bw is brighter 
than the DF. However at present it is not clear whether all GRBs are associated with the 
SNe or not (e.g., Price et al. 2002). In any cases, deep searches following the XRF will give 
us valuable information. 

If the DF associated with an XRF is observed, we will be able to estimate the Lorentz 
factor and the viewing angle of the jets. Let the typical frequency or the break energy of the 
DF (XRF) be v DF = 5 df i>q (vxrf = Sxrf^'o), where 5 = 1/7(1 — (3cos8 v ) is the Doppler 
factor. When ^ <C 1, 7 > 1 and (7#t,) 2 1, we can derive 5df ~ 1/(27) an d $xrf ~ 
27/ (l&v ) 2 - Since we assume that the XRF is the GRB observed from the off- axis viewing 
angle, we may use the observational consequence for the break energy 5GRBhu' ~ 200 £ keV, 
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where £ ~ 0.5-2 (Preece et al. 2000). In our model, Sqrb becomes ~ 27. Then, we obtain 
7 ~ 100^ 1 / 2 (/iz/ £ >ir/5eV)~ 1 / 2 . On the other hand, we can derive vdf/vxrf ~ (^/2) 2 , which 
implies that we can also estimate the viewing angle. 
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Fig. 1. — The light curves of the X-ray flash as a function of the normalized observed time 
T/(r /cP'j 2 ), where we adopt r /j3cy 2 = 1 sec. We choose jA9 = 5, «b = —1, (5b = —3, 
7*/q = 200 keV and D = 1 Gpc. The flux is proportional to D~ 2 . 
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Fig. 2. — The light curves of the delayed flash as a function of the normalized observed time 
T/(r /c/97 2 ), where we adopt r /[3cy 2 = 1 sec. We choose 7 = 100, jA9 = 5, = — 1, 
/9s = _ 3, 7z/q = 200 keV and D — 1 Gpc. Our jet model predicts that the flux of the delayed 
flash is almost constant (F ~ 2 x 10~ 19 -DQp C ergs s _1 cm -2 ) with both the observed time and 
the viewing angle. 
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Fig. 3. — The light curves of the X-ray flash afterglow in the UV band are shown by 
varying the viewing angle 9 V . We fix parameters as A6 = 0.05, n = 1cm -3 , p = 2.25, 
E = 2 x 10 54 ergs, e e = 0.1, Eb = 0.01, and D — 1 Gpc. Boxes represent the light curves of 
the delayed flash in the same band. We choose a canonical set of parameters as 7 = 100 and 
r /(3ey 2 = 1 sec. The light curve of the delayed flash does not depend on the viewing angle 
6 V so much. For comparison, we show the light curves of the delayed flash with one of the 
fiducial parameters changed. Note that all the flux is proportional to D~ 2 , and the flux and 
the duration of the delayed flash are proportional to (r //3c7 2 ) -1 and r /(3ey 2 , respectively. 
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Fig. 4. — The same as Fig. 3 but for A0 = 0.1 and E = 5 x 10 53 ergs. Note that the geometry 
corrected total energy (A9) 2 E/2 is not altered. 



